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We apply a contour deformation technique in momentum space to the newly developed Gamow 
shell model, and study the drip-line nuclei 5 He, 6 He and 7 He. A major problem in Gamow shell- 
model studies of nuclear many-body systems is the increasing dimensionality of many-body config- 
urations due to the large number of resonant and complex continuum states necessary to reproduce 
bound and resonant state energies. We address this problem using two different effective operator 
approaches generalized to the complex momentum plane. These are the Lee-Suzuki similarity trans- 
formation method for complex interactions and the multi-reference perturbation theory method. The 
combination of these two approaches results in a large truncation of the relevant configurations com- 
pared with direct diagonalization. This offers interesting perspectives for studies of weakly bound 
systems. 



I. INTRODUCTION 

We expect that present and proposed nuclear structure 
research facilities for radioactive beams will open new 
territory into regions of heavier nuclei. Such systems 
pose significant challenges to existing nuclear structure 
models since many of these nuclei will be unstable and 
short-lived. How to deal with weakly bound systems and 
coupling to resonant states is an open and interesting 
problem in nuclear spectroscopy. Weakly bound systems 
cannot be properly described within a standard shell- 
model approach since even bound states exhibit a strong 
coupling with the continuum. 

It is therefore important to investigate theoretical 
methods that will allow for a description of systems in- 
volved in such element production. Ideally we would like 
to start from an ab initio approach with the free nucleon- 
nucleon interaction and eventually also three-body inter- 
actions as the basic building blocks for the derivation of 
an effective shell-model interaction. The newly developed 
Gamow shell model offers such a possibility, see for ex- 
ample Refs. [11111111111111113. Similarly, the recent 
work on the continuum shell-model by Volya and Zelevin- 
sky |10| conveys similar interesting prospectives. Here 
we focus on the Gamow shell model, which has proved 
to be a powerful tool in describing and understanding 
the formation of multi-particle resonances within a shell- 
model formulation. Representing the shell-model equa- 
tions using a Berggren basis 0, El 0,0, 0,0, allows 
for a simple interpretation of multi-particle resonances 
in terms of single-particle resonances, as opposed to the 
traditional harmonic oscillator representation, where res- 
onances never appear explicitly. 



Although the Gamow shell-model approach is a pow- 
erful tool in this respect, there are major computational 
and theoretical challenges that need to be overcome if 
we aim at a realistic description of weakly bound and 
unbound nuclei. One of the challenges regarding the 
Gamow shell model discussed in Refs. 0, was the 
problem of choosing a contour in the complex A-plane 
that in the many-particle case selects the physical inter- 
esting states from the dense distribution of continuum 
states. In Refs. 0, Q the authors employ a "square- 
well" contour, which in the two-particle case separates 
the resonances from the complex-continuum states. In 
the case where more than two particles are present in 
the shell-model space, the resonant states mix with the 
complex continuum states, and an identification of the 
multi-particle resonances becomes difficult. 

In this work we consider as a test case the light 
drip-line nuclei 5,6 ' 7 He, and the formation of resonances 
in these nuclei starting from a single-particle picture. 
These nuclei have also been studied with a number of 
other methods, see for example Ref. |J_7| . and references 
therein. We construct a single-particle basis using the 
contour deformation method in momentum space, dis- 
cussed in detail in Ref. [lij , see also Ref. 0] for further 
references on complex scaling. We show that choosing a 
rotated plus translated contour in the complex plane, a 
large portion of the many-particle energy surface is free 
from complex continuum states. This choice of contour 
isolates the physical resonances, and allows for a clear 
distinction of many-particle resonances from the dense 
distribution of complex continuum states, also in the case 
when the number of particles exceeds two. 

The most severe problem and future challenge is that 
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the shell-model dimension increases dramatically for n > 
2 particles moving in a large valence space, this is what 
we henceforth refer to as the dimensionality problem. 
Using a technique such as the traditional Lanczos itera- 
tion method [l9j fails in Gamow shell-model calculations. 
Dealing with large real symmetric matrices, the Lanczos 
scheme is a powerful method when one wishes to calcu- 
late the states lowest in energy. In Gamow shell-model 
calculations there may be a large number of complex con- 
tinuum states lying below the physical resonances in real 
energy. In addition it is difficult to predict where the 
multi-particle resonances will appear after diagonaliza- 
tion. In Refs. 0, 0] this problem was circumvented by 
choosing a small number of complex continuum states in 
the single-particle basis. It was also pointed out that the 
results obtained were not converged with respect to the 
number of single-particle continuum orbits. In Ref. 
another approach was considered, where at most two 
particles where allowed to move in complex continuum 
states. This was based on the assumption that these 
configurations play the dominant role in the formation of 
many-particle resonances, and configurations where more 
than two particles move in continuum states could be ne- 
glected. 

Our aim in this work is to propose an effective in- 
teraction scheme which allows for a much larger num- 
ber of complex continuum states in the calculations, and 
in addition takes into account the mixing of configura- 
tions where all particles may move in complex contin- 
uum states. We show that if one aims at accurate cal- 
culations of the multi-particle resonances, the effect of 
all particles moving in the continuum may not always be 
neglected. Our choice of contour allows for a perturba- 
tive treatment of the many-particle resonances, and we 
propose a perturbation theory based scheme which com- 
bines the Lee-Suzuki similarity transformation method 
[2fl I2H I22I |23| and the so-called multi-reference pertur- 
bation method 0, HE US to account for couplings with 
configurations where all single-particles move in complex 
continuum states. 

Presently, Gamow shell-model calculations have been 
performed with phenomenological nucleon-nucleon in- 
teractions. A major challenge is to construct effective 
nucleon-nucleon interactions for drip-line nuclei starting 
from a realistic nucleon-nucleon interaction. In this pa- 
per we focus on the choice of contour and the dimen- 
sionality problem. The effective nucleon-nucleon inter- 
action adopted is purely phenomenological. However, 
the scheme we present, although implemented with a 
phenomenological nucleon-nucleon interaction, allows to 
define effective interactions computed with the complex 
scaled single-particle basis. The problem of construct- 
ing an effective interaction based on present interaction 
models for the nucleon-nucleon force will be considered 
in a forthcoming work. 

The outline of this work is as follows. Sec. ILTl gives a 
brief description of the contour deformation method in 
momentum space, and presents calculations of the en- 



ergy spectrum of the nuclei 5,6,7 He. Sec. IIIII presents 
first the Lee-Suzuki transformation method generalized 
to complex interactions. Thereafter we apply the similar- 
ity transformation method to the unbound nucleus 7 He 
and give a convergence study of its J 7 ' = 3/2j~ ground 
state resonance. Sec . IIVI gives a brief outline of the multi- 
reference perturbation method, and its application to the 
ground state of 7 He. In Sec. [S we present an effective 
interaction scheme, which combines the Lee-Suzuki sim- 
ilarity transformation and the multi-reference perturba- 
tion method, for calculation of multi-particle resonances 
in weakly bound nuclei. Sec. I VII gives the conclusions of 
the present study and future perspectives and challenges 
for Gamow shell- model calculations. 



II. THE GAMOW SHELL MODEL 

The newly developed Gamow shell model has proved 
to be a powerful tool in describing and understanding 
multi-particle resonances appearing in nuclei near the 
drip- lines. Here we discuss how two- and three-particle 
resonances are formed in 6 He and He, and how they are 
to be understood in terms of a single-particle picture. 
The specific choice of contour in the complex fc-plane 
makes it easy to identify and interpret the multi-particle 
resonances. In this section no truncations are made, all 
possible configurations within a model space are used in 
the shell model calculations. 



A. Berggren Basis in the Momentum 
Representation 

In Ref. 

we studied the contour deformation method 
applied to the momentum space Schrodinger equation. 
It was discussed and shown how the specific choices of 
contours based on the analytic structure of the poten- 
tial may allow for a unified description of bound, anti- 
bound (virtual) and resonant states. We will apply this 
method to obtain a single-particle Berggren basis for use 
in Gamow shell-model calculations. Here we briefly out- 
line the contour deformation method, and refer the reader 
to Ref. [16] for a more rigorous discussion. 

The analytically continued Schrodinger equation on a 
general inversion symmetric contour takes the form 

ft 2 2 f 

— k 2 il) n i{k) + - dqq 2 Vi(k,q)ip n i(q) = E n ii> n i(k). 

l[l 7T J c + 

(1) 

Here both k and q are defined on an inversion symmetric 
contour C + in the lower half complex fc-plane, resulting 
in a closed integral equation. The cigenfunctions consti- 
tute a complete bi-orthogonal set, normalized according 
to the Berggren metric [Tl|, [l^, El EH, Ell , namely 

1 = E \mma + [ dkk 2 \Mk))(^(k)\- (2) 
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In this work we construct a single-particle Berggren en- 
semble on a rotated plus translated contour, Cr+t, m 
the complex /c-plane, studied in detail in Ref. ^3\- The 
contour C^ +T is part of the inversion symmetric contour 
Cr+t — Cfl+T + C# +T displayed in Fig. ^ The com- 



with 




FIG. 1: Contour C+ 



L\ + L2 + L3 is given by the solid 



line, while the contour C R+T is given by the dashed line. The 
contour Cr+t = C R+T + C# +T is inversion symmetric. The 
single-particle spectrum which is exposed by this contour is 
marked by filled circles • and the excluded spectrum by open 
circles o. The full spectrum includes bound states (B), anti- 
bound (A), decay (D) and capture (C) resonant states. 

plete set of single-particle orbits defined by this contour 
will then include anti-bound, bound and resonant states. 
This basis serves as our starting point for Gamow shell- 
model calculations. 



v*Jk,k') = 



: exp(- 



4af 



k W f kk ' 

(4) 

where the subscripts Ij refer to the single-particle orbital 
and angular momentum quantum numbers I and j, re- 
spectively. The term Ii+i/2(z) is a Bessel function of the 
first kind with complex arguments. Fitting this poten- 
tial to reproduce the 5 He single-particle spectrum and 
phase-shifts results in g c = 47.4 MeV, g„i = 5.86 MeV 
and a c — a a i = 2.3 fm _1 . 

In the complex /c-plane the Gaussian potential diverges 
exponentially for |Im[fc]| > |Re[/c]|. If we apply the com- 
plex scaling technique, which consists of solving the mo- 
mentum space Schrodinger equation on a purely rotated 
contour, we get the restriction 6 < 7r/4 on the rotation 
angle. Even for smaller angles we may get a poor con- 
vergence, since the Gaussian potential oscillates strongly 
along the rotated contours. On the other hand, choosing 
a contour of the type C^ +T solves this problem, allowing 
for a continuation in the third quadrant of the complex 
/c-plane. Furthermore, it yields a faster and smoother de- 
cay of the Gaussian potential along the chosen contour. 

Since 5 He has only resonances in its spectrum, viz., no 
anti-bound states, there is no need for an analytic con- 
tinuation in the third quadrant of the complex /c-plane, 
as done in Ref. 16] for the free nuclcon-nuclcon interac- 
tion. We choose a contour of the type C^ +T rotated with 
6 = 7r/4 and translated with |Im[fc]| = 0.4sin(7r/4) « 
0.28 fm _1 in the fourth quadrant of the complex /c-plane. 
Figs. [21 and |3 give plots of the single-particle spectrum in 
5 He for the spin-orbit partners p 3 / 2 andp]y 2 respectively. 
We have used 50 integration points along the rotated Cr, 
and the translated GY parts of the contour C^ +T in the 
complex fc-plane. 



B. Single-Particle Spectrum of 5 He 

We consider first the unbound nucleus 5 He. This nu- 
cleus may be modeled by an inert 4 He core with a neu- 
tron moving mainly in the resonant spin-orbit partners 
P3/2 and pi/2- The J 71 ' = 3/2]" resonance, to be asso- 
ciated with the single-particle orbit P3/2, is experimen- 
tally known to have a width of T as 0.60 MeV while the 
J v = 1/2 j" resonance, associated with the single-particle 
orbit P1/2, has a large width r«4 MeV. For more infor- 
mation on these systems, see for example the recent re- 
view by Jonson [lj| • The core-neutron interaction in 5 He 
may be phenomenologically modeled by the SBB (Sack, 
Biedenharn and Breit) potential [23|. The SBB poten- 
tial is of Gaussian type with a spin-orbit term, fitted to 
reproduce the neutron - 4 He scattering phase shifts. In 
momentum space the SBB potential, which consist of a 
central part c and a spin-orbit term al, reads 



T 



, resonance in He 
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V v (k, k 1 ) = V t j(k, k') + (a-l )Vf 3 L {k, k') 



(3) 



FIG. 2: Plot of the P3/2 single-particle spectrum in 5 He for 
a Gaussian single-particle potential. The resonance is well 
located. The remaining points represent the non-resonant 
continuum. 

Table [I] gives the convergence of the p 3 / 2 and the 
P1/2 single-particle resonances as function of integration 
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FIG. 3: Plot of the pi/i single-particle spectrum in 5 He for 
a Gaussian single-particle potential. The resonance is well 
located. The remaining points represent the non-resonant 
continuum. 



points along the contour Cr+t- We observe that with 
12 points along the rotated path and 12 points along 
the translated line, one has a reasonable convergence of 
the resonance energy, giving in total 48 single-particle 
states for the valence space consisting of the Ij orbits 
{P3/21P1/2} with their pertinent momenta k defined by 
the number of mesh points. It is clear that if several 
particles were to move in this space, the dimensional- 
ity would become enormous. It is therefore important, 
even at the single-particle stage, to optimize the distri- 
bution of continuum states, in order that the main fea- 
tures of the system are reproduced with a small num- 
ber of single-particle resonances and complex continuum 
states. Notice also that the calculated width of the 1/2" 



TABLE I: Convergence of p 3 / 2 and P\ii resonance energies in 
5 He as function of the number of integration points Nr along 
the rotated Cr and Nt along the translated part Ct of the 
contour. Energies are given in units of MeV. 







,r = 


3/2" 


r = 


1/2- 


Nr 


N T 


Re[E] 


Im[E] 


Re[E] 


Im[E] 


10 


10 


0.752321 


-0.329830 


2.148476 


-2.912522 


12 


12 


0.752495 


-0.327963 


2.152992 


-2.913609 


20 


20 


0.752476 


-0.328033 


2.154139 


-2.912148 


30 


30 


0.752476 


-0.328033 


2.154147 


-2.912162 


40 


40 


0.752476 


-0.328033 


2.154147 


-2.912162 



resonance is somewhat larger (w 6 MeV) than the exper- 
imental value (« 4 MeV), see Ref. [l7| . 



C. Two-Particle Resonances in He 

Here we present results for the resonant spectra of 6 He. 
We employ again a shell-model picture with 6 He modeled 
by an inert 4 He core and two valence neutrons moving 
in the Ij orbits {£>3/2,Pi/2}i ignoring the recoil from the 



core. The model space consists then of all momenta k de- 
fined by the set of mesh points along the various contours, 
pertinent to these two Ij orbits. Using the single-particle 
wave functions for 5 He of Subsec lHBl we can in turn con- 
struct an anti-symmetric two-body wave function based 
on these single-particle wave functions, viz., 



JM 



(1,2), 



(5) 



a<b 



where the indices a, b represent the various single-particle 
orbits. Here <I>^ & (1,2) is an anti-symmetric two-particle 
basis state in the j — j coupling scheme. The sum over 
single-particle orbits is limited by a < b since we deal 
with identical particles only. The expansion coefficients 
fulfill the completeness relation 



1 



b ) 



(6) 



i<b 



and the two-particle Berggren basis forms a complete set 
1 = ^|^(1,2))(^(1,2)|. (7) 



a<b 



Here (I>^(1,2)| is the complex conjugate of ($^(1, 2)|. 
As an effective two-neutron interaction Vij we use a phe- 
nomenological interaction of Gaussian type, separable in 
r-i,rj, given by 

Vtffa.rj) = y cxp(-a 2 (r l 2 +r J 2 ))^(y A (z).y A (j)). (8) 



Two model spaces were considered. The first case in- 
cludes only the P3/2 single-particle orbit for various val- 
ues of the momentum k to be defined below. The second 
model space includes also the P1/2 single-particle orbit 
and its relevant momenta. For both model spaces we 
fit the interaction strength to reproduce the + binding 
energy in 6 He. We have observed that the position of 
the 2+ resonance in 6 He depends on the range a of the 
Gaussian interaction, even though the + ground state 
does not change with a. Unfortunately it turns out that 
for larger values of a the energy fit is better, but the 
convergence as function of meshpoints is poorer. In our 
calculations we have chosen a value of a which is a com- 
promise between a small number of mesh points along 
the contour and a reasonable good fit of the resonant 
energy spectra. This demonstrates that the two-particle 
resonant spectrum depends on the radial shape of the in- 
teraction and suggests that we should rather deal with 
an effective interaction derived from realistic models for 
the nucleon-nucleon interaction. The parameters used 
in our calculations are Vq = —5.315 MeV for the model 
space involving only the (P3/2) states and Vq — —4.549 
MeV for a model space consisting of both single-particle 
quantum states P3/2 and p\/2- We use a = 4.8 fm _1 for 
both model spaces. 

Figs. 01 and [S] show the + and 2 + energy spectrum, 
respectively, for 6 He after a full diagonalization of the 
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two-particle shell-model equation. Here the model space 
is defined by the p 3 / 2 and p 1 / 2 single-particle orbits. This 
model space yields a bound + state as well as a resonant 
+ state. Moreover, we obtain two resonant 2 + states. 
Observe that the choice of contour (C^ +T ) separates all 
physical relevant states from the dense distribution of 
complex continuum states in the energy plane. By this 
choice of contour the identification of multi-particle res- 
onances is fairly easy, and one may study the resonant 
trajectories as the interaction strength is varied. 



> 




FIG. 4: Plot of 0^ bound- and 2 resonant state in He for 
a model space consisting of the P3/2 and pi/ 2 single-particle 
orbits. The bound and resonant states are well located. The 
remaining points represent the non-resonant continuum. 
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FIG. 5: Plot of 2 + resonances in 6 He for a model space consist- 
ing of the p 3 /2 and p\/ 2 single-particle orbits. Both resonant 
states are well located. The remaining points represent the 
non-resonant continuum. 

The stability of the + and 2 + results as function of the 
number of mesh points is demonstrated in Tables [H] and 
IIIII Limiting first the attention to a model space consist- 



ing only of the P3/2 orbit, we note that with Nr = 12 in- 
tegration points along the rotated path Cr and Nt = 12 
points along the translated line C't, convergence is sat- 
isfactory, even with a total of 300 two-particle states. 



TABLE II: Convergence of the 0^ bound state energy in 6 He 
in terms of the number integration points Nr and Nt along 
the rotated Cr and the translated part Ct of the contour, 
respectively. The number N2 P gives the dimension of the two- 
particle anti-symmetrized basis. Here only p 3 / 2 single-particle 
orbits are included. Energies are given in units of MeV. 



Nr Nt 


N 2p 


Re[E] 


Im[E] 


12 12 


300 


-0.980067 


-0.000759 


20 20 


820 


-0.979508 


0.000000 


25 25 


1275 


-0.979509 


0.000000 



TABLE III: Convergence of the 2+ resonant state energy in 
6 He as function of the number of integration points Nr and 
Nt along the rotated Cr and the translated part Ct of the 
contour, respectively. The number N 2p gives the dimension 
of the two-particle anti-symmetrized basis. Here only p 3 / 2 
single-particle orbits are included, and energies are given in 
units of MeV. 



N R N t 


N 2p 


Re[E] 


Im[E] 


12 12 


300 


1.215956 


-0.267521 


20 20 


820 


1.216495 


-0.267745 


25 25 


1275 


1.216496 


-0.267745 



Tables El El and EH repeat the above convergence anal- 
ysis, but now employing a model space consisting of the 
P3/2 an d P1/2 single-particle orbits, and including also 
the results for the lowest-lying 6 He state with quantum 
numbers J w = 1 + . Increasing the model space brings 
several new features. We note in Table llVl that the first 
excited Oj" state is a resonance. However, the stabil- 
ity of the results as functions of the number of mesh 
points is comparable to that seen in Tables |H] and IIIII 
With approximately 12 mesh points we obtain results 
close to the converged ones. Similar conclusions apply 



TABLE IV: Convergence of Of bound and the 0]j~ resonant 
state energy in 6 He as function of the number of integra- 
tion points Nr and Nt along the rotated Cr and the trans- 
lated part Ct of the contour, respectively. The number N 2p 
gives the dimension of the two-particle anti-symmetrized ba- 
sis. Here the p 3 / 2 and Pi/ 2 single-particle orbits are included. 
Energies in units of MeV. 







r = 


= 0+ 


r -- 


= 0+ 


Nr Nt 


N 2p 


Re[E] 


Im[E] 


Re[E] 


Im[E] 


12 12 


600 


-0.980111 


-0.000497 


4.289194 


-3.882119 


20 20 


1640 


-0.979148 


-0.000000 


4.286186 


-3.882878 


25 25 


2550 


-0.979148 


0.000000 


4.286181 


-3.882876 



to the if resonance and the two lowest-lying 2 + reso- 
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nant states, see Tables IVl and fvTI for more details. We 

TABLE V: Convergence of the if resonance as function of the 
number of integration points Nr and Nt along the rotated 
Cr and the translated part Ct of the contour, respectively. 
The number N 2p gives the dimension of the two-particle anti- 
symmetrized basis. Here the p%/ 2 and P\/ 2 single-particle or- 
bits are included. Energies in units of MeV. 







r = 


= 1+ 


Nr N t 


N 2p 


Re[E] 


Im[E] 


12 12 


1128 


1.945539 


-2.920286 


20 20 


3160 


1.940263 


-2.930619 


25 25 


4950 


1.940266 


-2.930608 



TABLE VI: Convergence of the 2J and 2~f resonance energy 
in 6 He as function of the number of integration points Nr and 
Nt along the rotated Cr and the translated part Ct of the 
contour, respectively. The number N 2p gives the dimension 
of the two-particle anti-symmetrized basis. Here the ps/ 2 and 
P1/2 single-particle orbits are included. Energies in units of 
MeV. 







r = 


= 2+ 


r ■■ 


- 2+ 


Nr N t 


N 2p 


Re[E] 


Im[E] 


Re[E] 


Im[E] 


12 12 


876 


1.149842 


-0.203052 


2.372295 


-2.122474 


20 20 


2420 


1.150527 


-0.203060 


2.372818 


-2.123253 


25 25 


3775 


1.150527 


-0.203060 


2.372817 


-2.123254 



note that the experimental value for the width of the 
first excited J n = 2f is V w 113 KeV and the energy 
is Re[E] 2 + = 1797 KeV. Our simplified nucleon-nucleon 
interaction gives a qualitative reproduction of the data. 
In a future work we plan to include a realistic nucleon- 
nucleon interaction for studies of such systems. 

We end this subsection by analyzing the squared 
amplitude of the single-particle configurations 
\RR),\RC),\CC) of the 0+, 1+ and 2+ bound- and 
resonant wave functions. The results are shown in 
Tables lWrVTIdlliXl lXlandlXll The reason for doing this 
analysis is due to the fact that our single-particle basis 
consists of resonant and continuum single-particle orbits. 
By performing such an analysis we can disentangle 
the contribution from for example the non-resonant 
continuum. In these tables, \RR) stands for both 
single-particle orbits being a resonant single-particle 
orbit, \RC) means that one single-particle orbit is a 
resonant single-particle orbit and the other a non- 
resonant continuum single-particle orbit, while for \CC) 
both single-particle orbits are from the non-resonant 
single-particle continuum. All the results show that the 
configurations where both single-particles are resonant 
orbits, have the largest amplitude in the two-body wave 
function. It is also seen, that the configurations where 
both particles are in complex continuum states have a 
small effect on the formation of two-particle resonances 
in 6 He. This is a useful result which we will exploit 



below when we define effective interactions for smaller 
spaces. 

TABLE VII: Expansion coefficients of the Of bound state in 
6 He. The p 3 / 2 and p\/ 2 single-particle orbits define the model 
space. See text for further discussions. 





(pI 


ft) 


(Pi 


/ 2 ) 




Re[C 2 ] 


Im[C 2 ] 


Re[C 2 ] 


Im[C 2 ] 


\RR) 


1.10488 


-0.83161 


0.22620 


-0.16120 


\RC) 


-0.06036 


0.88137 


-0.19842 


0.22423 


\CC) 


-0.09716 


-0.04974 


0.02486 


-0.06305 



TABLE VIII: Expansion coefficients of the 0^ resonance in 
6 He. The pz/ 2 and p\n single-particle orbits define the model 
space. See text for further discussions. 





(Pi 


/ 2 ) 


Cp?/ 9 ) 




Re[C 2 ] 


Im[C 2 ] 


Re[C 2 ] Im[C 2 ] 


\RR) 


-0.01136 


-0.08003 


0.90189 0.33029 


\RC) 


0.04282 


-0.03939 


0.05966 -0.24478 


\CC) 


0.00617 


0.00494 


0.00082 0.02896 



TABLE IX: Expansion coefficients of the 1 + resonance in 
6 He. The Pz/ 2 and P1/2 single-particle orbits define the model 
space. See text for further discussions. 





(P1/2P3/2) 


(pi 


n) 


(pU) 




Re[C 2 ] Im[C 2 ] 


Re[C 2 ] 


Im[C 2 ] 


Re[C 2 ] Im[C 2 ] 


\RR) 


0.71068 -0.03739 








\RC) 


0.00381 -0.06948 


0.00016 


0.00003 


0.02224 0.01171 


\CR) 


0.20647 0.05208 


0.00037 


0.00071 


0.07067 0.03807 


\CC) 


-0.00984 -0.00149 


-0.00031 


-0.00009 


-0.00424 0.00585 



D. Three-Particle Resonances in 7 He 

Finally we consider the unbound nucleus 7 He, whose 
ground state ( J™ = 3/2~) is located sa 0.5 MeV above the 
6 He ground state, with a measured width T w 160 keV. 
Other continuum structures, with tentative spin assign- 
ments — 1/2 - , and J* = 5/2~, have been observed, 
see for example Ref. for an extensive review of the 
experimental situation. In this subsection we limit the at- 
tention to a model defined by the P3/2 single-particle or- 
bits only. Thus our single-particle basis from 5 He implies 
that only a J T = 3/2 resonance may be formed. The 
reason we do not include the pi / 2 single-particle orbits is 
that we aim at a diagonalization in the full space, tak- 
ing into account all complex continuum couplings. This 
model calculation will serve as a later reference. In the 
case of 24 mesh points in momentum space for the p 3 / 2 
single-particle quantum numbers Ij, the total dimension 
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TABLE X: Expansion coefficients of the 2^ resonance in 
6 He. The pz/2 an d P1/2 single-particle orbits define the model 
space. See text for further discussions. 



IP1/2P3/2J 



{pj/2 



Re[C 2 ] Im[C 2 ] Re[C 2 ] Im[C 2 ] 
\RR) 0.11394 -0.00494 0.96962 0.05539 
\RC) -0.00474 0.02531 -0.00178 -0.00018 
\CR) -0.02776 -0.03796 -0.05069 -0.02708 
\CC) 0.00229 -0.00772 -0.00089 -0.00282 



TABLE XI: Expansion coefficients of the 2^ resonance in 
6 He. The pz/2 an d P1/2 single-particle orbits define the model 
space. See text for further discussions. 



{P1/2P3/2) 



(P3/2) 



Re[C 2 ] Im[C 2 ] Re[C 2 ] Im[C 2 ] 
\RR) 0.88847 -0.03742 0.08911 -0.03742 
\RC) -0.04888 0.05674 -0.00104 -0.00055 
\CR) 0.06058 -0.03336 -0.00072 -0.02469 
\CC) 0.01131 -0.00447 0.00115 -0.00220 



d of the (J 7r = 3/2~) three-particle problem is d — 9224, 
if, in addition, we were to include 24 single-particle mo- 
menta for the P1/2 single-particle quantum numbers Ij, 
we would have roughly d ~ 40000 three-body configura- 
tions. 

In Refs. and the dimensionality problem was 
circumvented by choosing a small number complex con- 
tinuum states, typically of the order of five, although it 
was found that a larger number continuum states had to 
be included to obtain converged results. Our aim in this 
subsection is to study the full effect of a coupling to the 
complex continuum, and show that if one is to obtain 
accurate results, the effect of all particles moving in the 
continuum on the formation of multi-particle resonances 
may be important. In addition we include a larger num- 
ber of continuum states, in order to achieve a satisfactory 
convergence. As for 6 He, we construct a three-body wave 
function using the single-particle wave functions defined 
in 5 He. The three-body wave function is expanded in a 
three-particle anti-symmetric Berggren basis 

*i M (1, 2,3)= J2 C( a M b) y (a M b)c (l, 2, 3), (9) 

a<b<c 

where the completeness relation reads 

1= E l^) c (l,2,3))($f Q % c (l,2,3)|, (10) 



a<b<c 



with 



1 - ( C (<^b)c)Y 



(11) 



i<b<c 



The two-body nucleon-nucleon interaction is the same 
as that used for 6 He. Fig. Ogives the energy spectrum af- 
ter a full diagonalization of the three-particle shell- model 



equation. It is seen that the choice of contour in cal- 
culating the single-particle spectrum is again optimal if 
we wish to consider the fact that all physical interest- 
ing states are well separated from the dense distribution 
of complex scattering states. The J T = 3/2~ resonance 
appears at the energy E 3 / 2 - = —(0.12 + 0.12z) MeV. 

The plotted energy spectrum shows that the + and 2 + 
states in 6 He, and the 3/2~ state in 5 He, form complex 
thresholds in the energy spectrum of the J 77 = 3/2~ spec- 
trum in 7 He, see Tables . Ill ITT1 and Hill The physical inter- 
pretation of these three-particle states is, in the case of 
the 6 He thresholds, that two of the neutrons form either 
the + ground state or the 2 + resonant state, while the 
third neutron is moving in a complex continuum state. A 
diagonalization within the reduced space, where at most 
two particles move in continuum states gives the reso- 
nance energy — (0.14 + 0.16z) MeV, which shows that the 
effect coming from all particles moving in the continuum 
is not negligible, but small. 
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FIG. 6: Plot of the 3/2 complex energy spectrum of 7 He for a 
model space consisting of p 3 / 2 single-particle orbits only. The 



J w = 3/2 resonance is located at E- 
0.12221H) MeV. 



3/2- 



(0.120731 + 



Table IXIII gives the squared amplitudes of the vari- 
ous single-particle configurations in the 7 He ground state, 
{\RRR), \RRC), \RCC), \CCC)} , where R again labels a 
single-particle resonance and C a complex single-particle 
continuum orbit. It is seen that the most important con- 
figuration, as in the case of 6 He, is the one where all 
single-particles are in the p 3 / 2 single-particle resonant or- 
bit. The effect of configurations where all particles are 
in continuum states is small, which suggest that the cou- 
pling to configurations \CCC) may be taken into account 
perturbatively. This a feature we will exploit in Sees. IIIII 
andllVl 

In Fig. El we note that the J 71 " = 3/2~ ground state in 
7 He appears at energy of approximately 0.86 MeV above 
the ground state in 6 He, while the experimental value is 
at approximately 0.5 MeV. This discrepancy with exper- 
iment can be understood in terms of the configuration 
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\RRR), and the choice of interaction. Focusing on the 
first aspect and using coefficients of fractional parentage, 
we can rewrite the \RRR) configuration as 



need to extend its use to Gamow shell-model calcula- 
tions, implying a generalization to complex interactions. 
To achieve the latter we introduce first the two-body 
Schrodinger equation 



\ 1(^3/2)^3/2;^ = 3/2->-y|]0»3/a)2P3/a; J* = 3/2%+ V 12 )\^ n ) = £„|Vv, 



\(p 3 /2) 3 ;r = 3/2") 

From the geometry one may conclude that the ground 
state of 7 He bears much more resemblance with the 2^~ 
resonance than with the Of ground state of 6 He. In 
our calculations the 2+ resonance comes at an energy 
« (1.2-0.26z) McV, which is roughly 2.2 MeV above the 
Of ground state of 6 He, to be contrasted with the exper- 
imental value of w 1.8 MeV, see Fig. El This suggests 
that if we were to increase the attractive strength of the 
J w = 2 + interaction in 6 He and get a better agreement 
with the experimental value, the J T = 3/2 _ resonant 
ground state of 7 He would get closer to the experimental 
results. 



TABLE XII: Expansion coefficients of the J w = 3/2" ground 
state in 7 He. Here only p 3 / 2 single-particle orbits are included. 



(13) 



|pf/a) 




Re[C 2 ] 


Im[C 2 ] 


\RRR) 


1.295549 


-0.986836 


\RRC) 


-0.184544 


1.099729 


\RCC) 


-0.115738 


-0.110375 


\CCC) 


0.004733 


-0.002518 



III. 



EFFECTIVE INTERACTIONS FOR THE 
GAMOW SHELL MODEL 

The Lee-Suzuki Similarity Transformation for 
Complex Interactions 



The previous section served to introduce and motivate 
the application of complex scaling in studies of weakly 
bound nuclear systems. However, employing such a mo- 
mentum space basis soon exceeds feasible dimensionali- 
ties in shell-model studies. To circumvent this problem 
and to be able to define effective interactions of prac- 
tical use in shell- model calculations, we introduce effec- 
tive two-body interactions based on similarity transfor- 
mation methods. These interactions are in turn employed 
in Gamow shell-model calculations. We base our ap- 
proach on the extensive works of Suzuki, Okamoto, Lee 
and collaborators, see for example Refs. pel I2H I22I |23| . 
This similarity transformation method has been widely 
used in the construction of a effective two- and three- 
body interactions for use in the No-Core shell-model ap- 
proach of Barrett, Navratil, Vary and collaborators, see 
for example Refs. [2a. l29t f30L l31| and references therein. 
However, since the similarity transformation method has 
previously only been considered for real interactions, we 



here H$ includes the single-particle part of the Hamil- 
tonian, kinetic energy and an eventual single-particle in- 
teraction. The term V12 is the residual two-body inter- 
action. We then expand the exact wave function tp n in 
the anti-symmetric two-particle basis, generated by the 
single-particle basis of Hq, which corresponds to the ba- 
sis from the 5 He calculations of Subscc. Ill Bl Thereafter 
we choose a suitable single-particle model space p and its 
complement space q. These single-particle spaces define 
in turn our two- and many-particle model spaces 



p = X>p><fipi, 

ap 

and the complement space 

Q = Y1 I^qX^qI' 



(14) 



(15) 



where P is defined by both single-particle orbits being in 
the p— space, and the complement space Q is given by all 
two particle states where at least one particle is in the q- 
space. The anti-symmetric two-particle basis follows the 
Berggren metric 

(a'\a) = (a'*\a) = 8 a / <a , (16) 

and the projection operators fulfill the relations 

P 2 = P, Q 2 = Q, P T = P, (17) 

and 

Q T = Q, P + Q = l, PQ = 0. (18) 

We wish to construct an effective two-body interaction 
within the model space, reproducing in the P-space ex- 
actly the Np model space eigenvalues of the full Hamil- 
tonian. This can be accomplished by a similarity trans- 
formation 



H = e-^He" 



(19) 



where oj is defined by oj — QojP. It follows that 
oj 2 = oj 3 = ... = and e w = P + Q + oj. The two-body 
Schrodinger equation can then be rewritten in terms of 
a 2 x 2 block structure 



PHP PHQ\ 
QHP QHQ J 




(20) 



If PHP is to be the two-particle effective interaction, the 
decoupling condition PHQ = must be fulfilled. One 
may show that the decoupling condition becomes 



QHP + QHQoj - ujPHP - ojPHQuj = 0, 



(21) 
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with u> acting as a transformation from the model space 
P to its complement Q, viz., 



(ocQ\ip n ) = ^2(aQ\uj\ap)(a P \^ 71 



(22) 



There is however no unique solution for to. The effec- 
tive interaction generated in the model space depends 
on the Np exact solutions entering Eq. I|22|) . This is 
why the effective interaction generated by the similarity 
transformation method is often referred to as a state de- 
pendent effective interaction. The solution for lu may be 
obtained as long as the matrix (ap\ip n ) is invertible and 
non-singular. Based on this we choose those Np exact 
solutions with the largest overlap with the two-particle 
model space. With the solution u>, the non-symmetric 
effective interaction R is given by 



R = PHP - PH Q P = PV 12 P + PV 12 Quj. 



(23) 



ft would be preferable to obtain a complex symmetric 
effective interaction, in order to take advantage of the 
anti-symmctrization of the two-particle basis. This may 
be accomplished by a complex orthogonal transformation 

V cff = U- 1 (H a + V 12 )U -H , (24) 

where U is complex orthogonal and defined by 

U = exp(-S), S = arctanh(u - lu t ), (25) 

and 

U T U = UU T = 1, U T 



u- 



(26) 



Such complex orthogonal transformations preserve the 
Berggren metric x T x of any vector x £ {C n }. This fea- 
ture allows us to define a complex and symmetric effective 
two-body interaction 

V cS = {P+uj T uj) 1/2 {PHP + PHQuj)(P+Lu T ujy 



1/2. 



-H . 

In determining V e f{, one has to find the square root of 
the matrix A — (P + uj t lu). In the case of A being real 
and positive definite the method based on eigenvector 
decomposition gives generally a stable solution. Using 
the eigenvector decomposition, with Z representing an 
orthogonal matrix, D being a diagonal matrix composed 
of the eigenvalues, using A = ZDZ T , Z T Z = ZZ T = 1, 
D = (D) 1 /2(£>)i/2 j and 



A = ( ZD 1 ' 2 Z T 



ZD 1/2 Z T 



(28) 



we can write the square root of a matrix A as 

A 1/2 = ZD 1/2 Z T . (29) 

For a complex matrix A the procedure based on eigen- 
vector decomposition is generally numerically unstable. 
An approach suitable for complex matrices is based on 
properties of the matrix sign function. It can be shown 
that the square root of the matrix is related to the matrix 



sign function, see Ref. [32| for more details. In the case of 
A being complex and having all eigenvalues in the open 
right half complex plane, iterations based on the matrix 
sign function are generally more stable 








a' 




sign ^ 


J 





)■ 



A 1 ' 2 
A- 1 / 2 



(30) 



One stable iteration scheme for the matrix sign was de- 
rived by Denman and Beavers (3j| , as a special case of a 
method for solving the algebraic Riccati equation 



Yn = A, Z a =L 



Y k+1 = 



z, 



k+ 



i = \{Z k +Y^), fc = 0,l,2,. 



(31) 
(32) 

(33) 



and provided A has no non-positive eigenvalues this iter- 
ation scheme exhibits a quadratic convergence rate with 



tt-AVa, Z^A- 1 ' 2 as k 



oo. 



(34) 



In our calculations, convergence is typically obtained af- 
ter a few iterations. 



B. Gamow Shell-Model Studies of 7 He with the 
Similarity Transformation Method 

Here we apply the Lee-Suzuki similarity transforma- 
tion method to the Gamow shell- model calculation of the 
ground state of 7 He. The first problem is to define an 
optimal single-particle model space, which subsequently 
defines the two-particle model space, where the effective 
interaction is constructed. In Sec. [H] it was shown that 
the J v = 3/2~ ground state of 7 He has the 2 + reso- 
nance in 6 He as an important two-body configuration, see 
Eq. Q12JI. Based on this result, a viable starting point is 
to study the single-particle strengths in the 2 + resonance 
wave function. To understand the nature of two-particle 
resonances and how they are formed in a shell-model 
framework, it is natural to study and analyze the single- 
particle strengths in the two-particle wave function, and 
how they are distributed among the single-particle reso- 
nances and the various complex continuum orbits, given 
on a specific contour in the complex /c-plane. The single- 
particle density operator is given by 



N 



£hfctf)><&0')l, n = J2' 



(35) 



where TV is the total number of particles, i = {h,ji} 
labels the single-particle quantum numbers and i repre- 
sents the single-particle orbit. In the case of 6 He with 
an inert 4 He core, N = 2. Finding the probability, rij, 
that either particle 1 or particle 2 is in the single-particle 
orbit i = {h,ji}, we calculate the matrix element of hi 
with the two-particle resonance wave function 
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/ 1 \ 1/2 

, = (^(1,2)1^1^(1,2))^^^^.^ x 



Fig. [7| gives the real, imaginary and absolute values of 
the single-particle strengths among the complex contin- 
uum orbits, in the 2 + resonance in 6 He. The strengths 
are plotted as a function of the absolute value of the 
complex continuum energy. Observe that the continuum 
states near the 2 + resonance in 6 He have the largest 
strengths. This may be understood as an interference 
effect between the single-particle resonances and the con- 
tinuum orbits located closest in energy (momentum) to 
the single-particle resonance. When defining a single- 
particle model space, we choose the single-particle reso- 
nant and complex continuum orbits with the largest abso- 
lute value of the single-particle strength. With this recipe 
we have a consistent way of defining a single-particle 
model space, which forms the basis for constructing an 
effective interaction in the two-particle model space. 

In Figs. |H1 and |51 we show the convergence of the real 
and imaginary part of the J* — 3/2 _ resonance in He, 
as function of an increasing single-particle model space. 
For comparison we plot the results for a diagonalization 
within the model space using the "bare" interaction. It 
is seen that results with the effective interaction con- 
structed with the similarity transformation method con- 
verges much faster than results obtained with the "bare" 
interaction. We see that a satisfactory convergence is ob- 
tained with 10—11 single-particle Berggren states in the 
single-particle model space p from 5 He, corresponding to 
« 700 — 800 three-particle states Np. Compared with 
the full dimension of the three-particle problem, 9224, 
we have effectively reduced the dimension to about 8% 
of the full space. This is a considerable benefit which may 
allow us to extend the Gamow shell model with a corn- 



Here D®Q is the diagonal part and H ®Q the off-diagonal 
part of H®Q . In this form, we see that H° defines the 
unperturbed part while H 1 gives the perturbations to 
H°. In the first step we construct the complex orthogonal 
matrix </> which diagonalizes H pp . Here the columns of 
the matrix span the reference space P, and is a more 
convenient basis for perturbation expansions. 



I 

plex scaled single-particle basis to heavier systems and 
realistic effective interactions. However, we can further 
improve upon this approach by considering perturbative 
techniques as well. That is the topic of the next Section. 



IV. THE MULTI-REFERENCE 
PERTURBATION METHOD 

The Moller-Plesset multi-reference perturbation 
method has recently been revived in quantum chemistry, 
see for example Refs. [3 HE H3- with an emphasis 
on scattering theory and electron decays in many-body 
systems. Here we only give a brief outline of the method, 
and refer the reader to Refs. [2~H I2M |26| for further 
details. 

The basic idea of the multi-reference perturbation 
method is to first diagonalize within a small space (refer- 
ence space) , and then add a perturbation to the reference 
states by taking into account excitations from the refer- 
ence space to the complement space. 

First we define a suitable ./V-particle reference (model) 
space P which describes most of the many-body correla- 
tions of the system, and hopefully gives a weak coupling 
with the complement space Q. The ./V-body problem may 
then be written as a block structure 

Thereafter we divide the full Hamiltonian in two parts 



(38) 

I 

sion in energy, and define M = tfi T H p Q which gives the 
orthogonal transformation of the coupling block with re- 
spect to the reference states cf>. Using intermediate nor- 
malization, the energy corrections up to third order for a 
given state <pi in the reference space, may then be shown 
to be, 



h pp h pq 



H pp 
Dfl« 



H p Q 
ffQP jjQQ 



H» + H 



Secondly, we perform a standard perturbation expan- 



Ef = $H pp <k, El = 0, 



(39) 
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FIG. 7: Plot of the P3/2 single-particle complex continuum strength m in the two-particle resonance wave function 2+ in He. 
The solid line gives the real part, the dotted line the imaginary and the dash-dotted line the absolute value of the strengths. 
The filled circles give the actual location of the complex continuum states in absolute value of energy. 
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FIG. 8: Convergence of the real part of the J w = 3/2 ~ reso- 
nance in 7 He for a space defined by the p 3 / 2 single-particle or- 
bits only. The abscissa represents the number of three-particle 
model space configurations N P while n sp represents the to- 
tal number of single-particle momenta for the p 3 / 2 single- 
particle quantum numbers Ij. The solid line represents the 
effective interaction generated by the Lee-Suzuki similarity 
transformation method, and the dashed line is obtained us- 
ing the the bare interaction and the same number of three- 
body configurations. The 3/2~ resonance is located at E 
= -(0.120731 + 0.12221H) MeV. The horizontal line is the 
real energy obtained in the full space of three-body configu- 
rations. 



FIG. 9: Convergence of imaginary part of J n = 3/2 ~ 
resonance in 7 He number of three-particle model space 
configurationsiVp while n sp represents the total number of 
single-particle momenta for the p 3 / 2 single-particle quantum 
numbers Ij. The solid line represents the effective inter- 
action generated by the Lee-Suzuki similarity transforma- 
tion method, and the dashed line is obtained using the the 
bare interaction and the same number of three-body con- 
figurations. The J"" = 3/2~ resonance is located at E 
= -(0.120731 + 0.122211i) MeV. The horizontal line is the 
imaginary energy obtained in the full space of three-body con- 
figurations. 



here Mi j is the dot product of the vector cfif with column 

n ^2 j of the coupling block H PQ . Observe that there is no 

Ef = — g — % ' 3 , (40) first order correction in the energy, meaning that it has 

j=N P +i 1 ~ i'j been accounted for by the reference states and energies. 

N M- H Q ®M T ^ n tne a PP uca ti° n of the multi-reference perturbation 

E 3 = — — — — i ^ k (411 method to the calculation of multi-particle resonances 

j.k=N P +i ^ » 3:3 " - i n k,k) in Gamow shell-model calculations, we have to define 
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a multi-particle model space which describes most of 
the many-body correlations. Based on our knowledge 
from the two-particle system, we assume a good choice 
for the complement space to consist of configurations 
where all particles move in complex continuum states. 
One would expect that the most important configura- 
tions of the multi-particle resonance, are configurations 
which include single-particle resonances. In our calcu- 
lation of 7 He, the three-particle model space is chosen 
by studying the squared amplitudes of the three-particle 
configurations given in Table IXTTl There we see that the 
amplitudes of configurations where all particles move in 
continuum states are small. In Refs. [(| and 0, where the 
Helium isotopes 6_9 He were studied within the Gamow 
shell-model formulation, the authors reached similar con- 
clusions. Note also that we use the bare two-body inter- 
action of Eq. (JSJ . This is to be contrasted to the method 
outlined in the next section where we use the Lee-Suzuki 
transformation in order to define an effective interaction. 

The three-particle model space, and corresponding 
complement space, used in our calculations is defined by 



\RRR), \RRC), \RCC), 
r = I Re(e a + e b + e c ) < E cutl 
Im(e a + e b + e c ) > -E cui 



Q = l-P, (42) 
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FIG. 11: Convergence of the real part of the J™ = 3/2 _ 
resonance in 7 He, as the dimension of three-particle model 
space increases with increasing cutoff in energy. The cutoff in 
energy is increased in steps of 5 MeV, i.e. E cut — 0, 5, 30, 
and given by the filled circles. The horizontal line is the real 
part of the J n = 3/2~ resonance located at E = —(0.120731 + 
0.12221H) MeV. The dashed-dotted line is the zeroth order 
energy, the dashed line represents the second-order energy 
and the solid line is the third-order energy. 



where at most two particles move in continuum states. 
We have also introduced a rectangular cutoff in the com- 
plex energy plane, since we assume that three-particle 
configurations high in the energy play a minor role on 
the formation of low-lying resonances. Fig. 1101 gives a 
plot of the unperturbed three-particle spectrum where 
at most two particles move in complex continuum states, 
and three different cut-offs in energies and corresponding 
model spaces are shown. 
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FIG. 10: Three choices of the model space used in the multi- 
configuration perturbation method. The three-particle model 
space states are constructed such that at most two particles 
move in the non-resonant continuum. 



Figs. ^2 an d El snow the convergence of the real and 
imaginary part of the three-particle resonance energy in 
the multi-reference perturbation method up to third or- 
der. The model space used here is given in Eq. 1421 
and the calculations were done for the increasing cut- 
offs in energy, E cut — 0,5, ...,30 MeV. Here the conver- 
gence is plotted with respect to the number of three- 
particle model-space states Np for each energy cutoff. 
We see that a satisfactory convergence is obtained with 
Np w 400, corresponding to the energy cutoff E cut = 20 
MeV. As expected, we see that excitations of model space 
configurations located above E cut ss 5 MeV yield small 
contributions to the second- and third-order corrections 
to the resonance energy. Observe that the second- and 
third-order terms converge at the same number of model 
space states, which indicates that second-order correc- 
tions in energy are seemingly sufficient for our applica- 
tions. This is also an advantage from a numerical point 
of view. In second order one has to store only the di- 
agonal part of the block H®® , while in third order the 
complete block H®® has be stored, which may be ex- 
tremely large in many cases. The zeroth order energy, 
which corresponds to diagonalization within P, does not 
saturate at the exact resonance energy with increasing 
Np, which again shows that possible couplings with the 
Q— space have to be accounted for, if one aims at accu- 
rate calculations. 

Summing up these results, we see that we obtain stable 
results with approximately Np 400 three-body config- 
urations within the multi-reference perturbation method, 
while the similarity transformation method of Sec. Ifffl 
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FIG. 12: Convergence of the imaginary part of the J n = 
3/2 ~ resonance in 7 He, as the number of three-particle model 
space increases with increasing cutoff in energy. The cutoff in 
energy is increased in steps of 5 MeV, i.e. E cut = 0, 5, 30, 
and given by the filled circles. The horizontal line is the the 
imaginary part of the J w = 3/2~ resonance located at E 
= -(0.120731 + 0.12221H) MeV. The dashed-dotted line is 
the zeroth order energy, the dashed line represents the second- 
order energy while the solid line is the third-order energy. 



gives stable results for Np » 800 three-body configura- 
tions for the same problem. The question now is whether 
we can marry these two approaches in our quest for 
smaller Gamow shell-model spaces. This is the topic of 
the next section. 



V. EFFECTIVE INTERACTION SCHEME FOR 
THE GAMOW SHELL MODEL 

In the previous sections it was shown that the Lee- 
Suzuki similarity transformation and the multi-reference 
perturbation method may be used in the Gamow shell 
model in order to account for the most important cor- 
relations of for example a multi-particle resonance. Al- 
though the dimensionality of the problem derived either 
from the similarity transformation method or the multi- 
reference perturbation method was small compared to 
the full problem, the dimensionality may still be a severe 
problem when dealing with more than three particles in 
a big valence space. 

The drawback of the multi-reference perturbation 
method is that one has to store extremely large matrices 
jjQQ if one wishes to go beyond second order in pertur- 
bation theory. In the similarity transformation method 
one does not have to deal with H®®, as couplings with 
the Q-space states have been dealt with, in practical cal- 
culations at least at the two-body level. Going to systems 
with larger degrees of freedom, the P-space may never- 
theless, at the converged level, be too large for our brute 
force diagonalization approach. 



The aim of this section is to propose an effective inter- 
action and perturbation theory scheme for the Gamow 
shell model. This approach combines the similarity 
transformation method and the multi-reference perturba- 
tion method, so that hopefully multi-particle resonances 
where several particles move in large valence spaces, may 
be calculated without a diagonalization in the full space. 
Our algorithm is as follows 

1. Choose an optimal set of n sp single-particle orbits, 
which in turn define two-body Pi v and many-body 
spaces. In our case these single-particle orbits are 
defined by selected states in 5 He. 

2. Construct a two-particle effective interaction by 
the Lee-Suzuki similarity transformation method 
within the two-particle model space Pi v . Such di- 
agonalizations can be done for very large spaces, 
see for example Refs. jH |H H3, [3lj . 

3. The next step is to divide the multi-particle model 
space P in two smaller spaces P' and Q' , where 
P = P' + Q' and N P = N P , + N Q > . The choice 
of P' should be dictated by our knowledge of the 
physical system. As an example, one may con- 
sider those single-particle configurations within the 
P— space that play the dominant role in the forma- 
tion of the multi-particle resonance. The number 
Np = Npr + Nqi represents the total number of 
many-body configurations within the P— space. 

4. Now that we have divided the P-space in two sub- 
spaces P' and Q', we use for example the multi- 
reference perturbation method to account for ex- 
citations from the P'— space to the Q'— space 
to obtain energy corrections to a specific order. 
Increase the size of the P' — space until conver- 
gence is obtained. In the case Npr = Np and 
Nq' = Np — Npr = the multi-reference pertur- 
bation expansion terminates at zeroth order, and 
corresponds to a full diagonalization within the 
P— space. Another option is to use for example the 
coupled cluster method as exposed in Refs. |34ll35|. 

5. Start from top again with a larger set of single- 
particle orbits, and continue until a convergence 
criterion is reached. 

We illustrate these various choices of model spaces in the 
following two figures. Fig. ^| defines our model space 
for the Lee-Suzuki similarity transformation at the two- 
body level. This corresponds to steps one and two in the 
above algorithm. The set of single-particle orbits defines 
the last single-particle orbit in the model space n sp . Note 
that we could have chosen a model space defined by a cut 
in energy, as done by the No-Core collaboration, see for 
example Refs. [2^, |2jJ 13(1 El- These examples serve 
just to illustrate the algorithm. Fig. [21 demonstrates 
again a possible division of the space into the full model 
space Np and a smaller space Npr. Again, this figure 
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basis. 

We define five different model spaces P, given by the 
total number of three-body configurations Np. The num- 
ber of single-particle orbits and three-body states are 
listed in Table. IXIIII The single-particle model space, 
defining P, is constructed according to the prescription 
outlined in section 11111 The reference space P' which 
defines a proper subset of each model space P, is again 
defined by 



\RRR),\RRC),\RCC) 
P'i = { Re(e a + e b + e c ) < E cut 
Im(e a + e b + e c ) > --E C ut 



C Pi. 



(43) 



FIG. 13: Possible definition of the two-body exclusion 
operator Qi v = 1 — P^ p used to compute the Lee-Suzuki 
similarity transformation and its effective interaction at 
the two-body level. The border of the model space is 
defined by the last single-particle orbit n sp . 



N P 



The reader should note that for each space Pi, P2 and 



TABLE XIII: Five different P— spaces defined for increasing 
number of single-particle model space orbits n sp consisting of 
the Ij configuration p 3 / 2 . The number Np gives the dimension 
of the three-particle model space P for = 3/2~ with a full 
dimensionality with n 3p = 24 of Np = 9224. 



p 


Pi P2 


P 3 


Pi Ps 


TLsp 


8 10 


12 


14 16 


N P 


344 670 


1156 


1834 2736 



so forth listed in Table IXIHI we can compare the results 
from this perturbative analysis with those from the exact 
diagonalization done in these spaces. This is shown in 
Tables IXIVIXVI IXVll IXVlTl and IXVIIII 



N P , N P 



FIG. 14: Possible definition of many-body space Np 
and reduced space N' P . 



serves only the purpose of illustrating the method. In 
our actual calculations we define the smaller space Np> 
via an energy cut in the real and imaginary eigenvalues 
and selected many-body configurations. 

In summary, defining a set of single-particle orbits in 
order to construct the two-body and many-body model 
spaces, we obtain first an effective two-body interac- 
tion in th e sp ace p2 P by performing the Lee-Suzuki 
[20I I2H I22L |23| transformation. This interaction and the 
pertinent single-particle orbits are then used to define a 
large many-body space. It is therefore of interest to see if 
we can reduce this dimensionality through the definition 
of smaller spaces and perturbative corrections. 

We present here as a test case, the calculation of the 
J n = 3/2 ~ three-particle resonance within the pertur- 
bation scheme outlined above. 24 single-particle orbits 
for the Ij configuration p 3 / 2 are included, giving a total 
dimension of d — 9224 for the J — 3/2 three-particle 



TABLE XIV: Resonance energy to second (E 2 ) and third 
order (E 3 ) in the multi-reference perturbation expansion, for 
the model space Pi given in Table IXIIII The subspaces P'i 
are defined for different energy cutoffs, increased in steps of 
10 MeV. In the last line we give the exact energy within Pi. 
Energies are given in units of MeV. 



N Pl = 


344 


N P 


-max , 


113 


N pll iV Q / Pcut 


Re[E 2 ] 


Im[E 2 ] 


Re[E 3 


Im[E 3 ] 


1 343 
113 231 10 


0.066 
0.041 


0.322 
-0.075 


0.606 
0.041 


0.088 
-0.076 


Exact within Pi 


0.042 


-0.076 







As the number of reference states Np, increases with 
increasing cutoff in energy P cu t , one reaches a maximum 
of reference states A?p/ max within each P— space. From 
the definition of the reference space P' in Eq. (|43J) . it 
will never coincide with the P— space as one exhausts the 
number of configurations \RRR),\RRC),\RCC) within 
P, since one by definition never includes the configura- 
tions \CCC), i.e. {|P')} C {|P)}. The perturbation 
scheme for a reference space P' given by Eq. H43|) . will 
therefore only yield convergent results as long as our as- 
sumption that the configurations \CCC) play a minor 
role compared to the reference states. Although the con- 
figurations \CCC) turn out to play a minor role for the 
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TABLE XV: Resonance energy to second (E 2 ) and third 
order (E 3 ) in the multi-reference perturbation expansion, for 
the model space Pi given in Table IXTTTI The subspaces P' 2 
are defined for different energy cutoffs, increased in steps of 
10 MeV. In the last line we give the exact energy within P2. 
Energies are given in units of MeV. 





N P2 = 


670 


N P 




181 




E C ut 


Re[E 2 ] 


Im[E 2 ] 


Re[E 3 ] 


Im[E 3 ] 


1 669 





-0.053 


0.357 


0.562 


0.059 


157 513 


10 


-0.078 


-0.110 


-0.079 


-0.110 


181 489 


20 


-0.082 


-0.110 


-0.083 


-0.110 


Exact within Pi 


: -0.081 


-0.110 







TABLE XVI: Resonance energy to second (E 2 ) and third 
order (E 3 ) in the multi-reference perturbation expansion, for 
the model space P3 given in Table IXTTTI The subspaces P' 3 
are defined for different energy cutoffs, increased in steps of 
10 MeV. In the last line we give the exact energy within P3. 
Energies are given in units of MeV. 



Np 3 = 


1156 


N P 




265 


N P , N q , r 


E cu t 


Re[E 2 ] 


Im[E 2 ] 


Re[E 3 ] 


Im[E 3 ] 


1 1155 





-0.099 


0.378 


0.561 


0.050 


205 951 


10 


-0.114 


-0.134 


-0.114 


-0.133 


265 891 


20 


-0.117 


-0.130 


-0.118 


-0.130 


Exact within P3 


: -0.116 


-0.130 







states we have considered in this work, there is no a 
priori reason for this to be the case when considering 
other multi-particle resonances. If no convergence is ob- 
served, one should simply choose another reference space 
P' , based for example on the single-particle model space, 
see Figs. El and fill 

Figs. El an d El gives plots of the real and imaginary 
part of the resonance energy to third order in the multi- 
reference perturbation expansion for the different model 
spaces considered above. From the plot one concludes 
that convergence is obtained for a small number of refer- 
ence states N P , ~ 350 - 400. 



TABLE XVII: Resonance energy to second (E 2 ) and third 
order (E 3 ) in the multi-reference perturbation expansion, for 
the model space P4 given in Table IXTTTI The subspaces P\ 
are defined for different energy cutoffs, increased in steps of 
10 MeV. In the last line we give the exact energy within P4. 
Energies are given in units of MeV. 





N Pi = 


1834 


N P 


1 m aa . 


365 




N Q ' 4 


E cu t 


Re[E 2 ] 


Im[E 2 ] 


Re[E 3 ] 


Im[E 3 ] 


1 


1833 





-0.134 


0.397 


0.532 


0.026 


253 


1581 


10 


-0.155 


-0.160 


-0.130 


-0.141 


347 


1487 


20 


-0.119 


-0.127 


-0.120 


-0.126 


365 


1469 


30 


-0.122 


-0.123 


-0.123 


-0.124 


Exact within P4 


: -0.121 


-0.124 







TABLE XVIII: Resonance energy to second (E 2 ) and third 
order (E 3 ) in the multi-reference perturbation expansion, for 
the model space P5 given in Table 1X1111 The subspaces P'5 
are defined for different energy cutoffs, increased in steps of 
10 MeV. In the last line we give the exact energy within P5. 
Energies are given in units of MeV. 





Np 5 = 


2736 




419 






E cu t 


Re[E 2 ] 


Im[E 2 ] 


Re[E 3 ] 


Im[E 3 ] 


1 


2735 





-0.137 


0.399 


0.530 


0.024 


253 


2483 


10 


-0.159 


-0.160 


-0.131 


-0.141 


347 


2389 


20 


-0.120 


-0.129 


-0.118 


-0.125 


409 


2327 


30 


-0.122 


-0.122 


-0.122 


-0.122 


419 


2317 


40 


-0.122 


-0.122 


-0.123 


-0.122 


Exact within P5 


: -0.121 


-0.122 








FIG. 15: Convergence of the real part of the J 71 = 3/2~ reso- 
nance energy in 7 He within the perturbative scheme outlined 
in the text, for the different model spaces P given in Ta- 
ble IXfffl The open circles along the different solid lines gives 
the calculations within each Pi . Npi gives the number of ref- 
erence states in P' , which is a subspace of P. The horizontal 
line is the the real part of the J 71 = 3/2 ~ resonance located 
at E = -(0.120731 + 0.122211*) MeV. 



In the approach considered above, the dimension of the 
Q'— space is considerably smaller than the dimension of 
the complement space Q = 1 — P, which makes it much 
less time and memory consuming to compute the matrix 
elements of . We have seen from the above calcula- 

tions that a termination of the perturbation expansion at 
second order compares well with the rate of convergence 
for the third-order expansion. This makes it numerically 
feasible to treat systems where several particles move in a 
large valence space, within perturbative scheme outlined 
above. 

We conclude this work by applying our scheme to 
the calculation of the three-particle resonances in 7 He, 
where 24 single particle states for each of the Ij single- 
particle states pi/2 and p 3 / 2 are included. The Hamil- 
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Level Scheme for 5 7 He 



r = l2.3 




FIG. 16: Convergence of the imaginary part of the J 71 = 3/2~ 
resonance energy in 7 He within the perturbative scheme out- 
lined in the text, for the different model spaces P given in 
Table 1X1111 The open circles along the different solid lines 
gives the calculations within each Pi. N P i gives the num- 
ber of reference states in P', which is a subspace of P. The 
horizontal line is the the imaginary part of the J n = 3/2~ 
resonance located at E = -(0.120731 + 0.12221H) MeV. 



tonian for the J n = l/2 _ ,3/2 _ and J* = 5/2" 
states for 7 He has dimensions Np — 29648, 38896 
and Np = 27072, respectively. The main compo- 
nent of the resonant wave functions turns out to be 
the \RRR) configuration, which gives the allowed cou- 
plings and corresponding unperturbed energies and wave 
functions | (P3/2) 3 ; J* — 3/2~~) with energy E = 
(2.2560 - 0.9840i) MeV, \p 1/2 (p 3/2 )l ; J* = l/2~ ) with 
energy E Q = (3.6581 - 3.5681*) MeV, \p 1/2 (p 3/2 ) \ ; J* = 
3/2^, 5 / 2 i > with energy E = (3.6581 - 3.568H) 
MeV and | (pi/ 2 ) Q Ps/ 2 ;J n = 3/2") and energy E Q = 
(5.0602 — 6. 1522i). We report here only the converged 
results for the lowest-lying 7 He resonances. They are 
2?(3/2~) = (0.02-0.08z) MeV, 22(1/2") = (0.39-3.98i) 
MeV, 2?(3/2J) = (2.43-1.95i) MeV, E(5/2~) = (2.75- 
0.89i)MeV, and E(3/2^) = (3.85 - 3.06i) MeV. Using 
our combination of the Lee-Suzuki similarity transfor- 
mation and the multi-reference perturbation method, re- 
sults close to the exact ones where obtained with approx- 
imately Npi ~ 1400 — 2000 three-particle configurations. 
This is a considerable reduction compared with the full 
dimensionalities listed above. 

Fig. El displays the calculated energy levels for the 
nuclei 5_7 He within our model. The unperturbed energy 
levels for 6 He and 7 He are also shown, and serve to illus- 
trate how the two- and three-particle resonances develop 
when the nucleon-nucleon interaction from Eq. © is in- 
cluded. There are several interesting features which can 
be seen from Fig. El The + - and 2 + -states in 6 He are 
formed within our model due to a strong pairing effect 



FIG. 17: Energy levels for the nuclei 5,6,7 He for model space 
consisting of the Ij single-particle states P\n and P3/2- The 
unperturbed energy levels for 6 He and 7 He are shown to the 
left in each case. 



between the two neutrons moving in equivalent orbits. 
When the nucleon-nucleon interaction is included, we ob- 
serve that the + and 2 + unperturbed energies gain an 
additional attraction. On the other hand, the 1 + state 
in 6 He exhibits weak pairing effects, the main component 
of the 1 + wave function is still the configuration \RR), 
but here the dominant contribution comes from a cou- 
pling between the pi/ 2 and P3/2 single-particle resonances 
in He. Adding the nucleon-nucleon interaction, does 
not change the unperturbed energy level significantly, see 
Fig. El The calculations done within our crude model of 
7 He, may suggest that this unbound nucleus has an even 
richer continuum structure than proposed in the recent 
review of Jonson, see Ref 01 ■ ^ n R- e f- E3 two excited 
states with tentative spins J T = 1/2 and = 5/2~ 
where reported to exist above the ground state in 7 He. 
The level spacings relative to the 7 He ground state were 
reported to be 0.57 MeV and 2.87 MeV for the J n = 1 /2" 
and J T = 5/2 states, respectively. The main decay 
channel of the J T = 5/2~ resonance at 2.87 Mev is a+3n. 
From this decay channel, Jonson Ref. 17] concluded that 

the configuration \p\/ 2 (jp 3 / 2 ) 2 J n = 5/2~) is the most 
probable one. In our calculations we find a rich con- 
tinuum structure in the energy region Re [25] w 3 MeV 
above threshold. From Fig. ll7l it is seen that the resonant 
three-particle states J"" = 3/2jj~, 5/2jj~ and J v = 3/2jj~ are 
located rather close in real energy. Although the widths 
vary from 2 — 6 MeV, this observation raises the question 
of whether these structures may be observed, and which 
spin assignements and the nature of the experimentally 
observed structures around 2.87 MeV in 7 He are. Fur- 
ther it is seen that the real part of the J*' = l/2~ reso- 
nance changes strongly, and moves towards the threshold 
when the nucleon-nucleon interaction is included. How- 
ever, these results must be gauged with the fact that we 
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are using a purely phenomenological nucleon-nucleon in- 
teraction model. The inclusion of a realistic interaction 
is the topic for a future work. The main issue here was 
to demonstrate how to derive effective interactions for 
the Gamow shell model, with a considerable reduction in 
dimensionality. 

VI. CONCLUSION AND FUTURE 
PERSPECTIVES. 

In this work we have applied the contour deformation 
method in momentum space, with a single-particle ba- 
sis in momentum space serving as a starting point for 
Gamow shell-model calculations. The resonant spec- 
tra of the drip-line nuclei 5_7 He have been studied 
and described using phenomenologically derived nucleon- 
nucleon interactions. It was illustrated that the choice of 
contour gives a good convergence for various resonant 
states, and in addition allows for a clear distinction be- 
tween all physical states and the remaining complex con- 
tinuum states. The main purpose of this work was to 
propose an effective interaction scheme for Gamow shell- 
model calculations. One of the most severe difficulties re- 
garding Gamow shell-model calculations is the dramatic 
growth of the shell-model dimension when dealing with 
several valence particles moving in a large shell-model 
space. This dimensionality problem is even more severe 
than the harmonic oscillator representation used in tradi- 
tional shell-model equation studies, In the Berggren rep- 
resentation a large number of complex-continuum states 
has to be included as well. The clear distinction of 
the unperturbed resonances from the dense distribu- 
tion of complex continuum states, allows for a pertur- 
bation treatment, when configuration mixing is taken 



into account. For perturbation expansions to converge, 
the unperturbed states have to be well separated from 
the Q-space states, or else the propagators will contain 
poles which make a perturbative treatment difficult. It 
has been shown that the Lee-Suzuki similarity transfor- 
mation combined with the multi-reference perturbation 
method, reduces the full problem to about 3 — 4%. 

Treating the many-particle problem in some perturba- 
tion scheme, we need to define a reference (model) space 
which describes most of the many-body correlations. The 
method and scheme outlined here, allows for a perturba- 
tive treatment of many-body states in which anti-bound 
states may play an important role, such as in the drip- line 
nuclei 11 Li. 

As has been pointed out, the location of multi-particle 
resonances depends on the effective interaction used be- 
tween valence nucleons. The next step is to derive a 
realistic effective interaction for Gamow shell-model cal- 
culations, and self-consistent Hartree-Fock single-particle 
energies for loosely bound nuclei, starting from a realis- 
tic nucleon-nucleon force. Using the Berggren represen- 
tation may give an underlying understanding of many- 
body resonances from a microscopic point of view. More- 
over, in our algorithm of Sec. we employed the multi- 
reference perturbation method. Our future plans in- 
volve replacing this method by the Coupled Cluster ap- 
proaches, as discussed in Refs. j^Hll- 



Acknowledgments 

Support by the Research Council of Norway is greatly 
acknowledged. 



Nazarewicz, and M. Ploszajczak, 
M. Ploszajczak, and J. Ro- 



N. Michel, W 
nucl-th/0407110 
N. Michel, W. Nazarewicz, 
tureau, nucl-th/0401036 

J. Dobaczewski, N. Michel, W. Nazarewicz, 
M. Ploszajczak, and M. V. Stoitsov, |nucl-th/040 1(534] 
R. J. Liotta, E. Maglione, N. Sandulescu, and T. Vertse, 
Phys. Lett. B367, 1 (1996). 

R. Id Betan, R. J. Liotta, N. Sandulescu, and T. Vertse, 
Phys. Rev. C 67, 014322 (2003). 

N. Michel, W. Nazarewicz, M. Ploszajczak, and K. Ben- 

naceur, Phys. Rev. Lett. 89, 042502 (2002). 

N. Michel, W. Nazarewicz, M. Ploszajczak, and 

J. Okolowicz, Phys. Rev. C 67, 054311 (2003). 

R. Id Betan, R. J. Liotta, N. Sandulescu, and T. Vertse, 

Phys. Rev. Lett. 89, 042501 (2002). 

R. Id Betan, R. J. Liotta, N. Sandulescu, and T. Vertse, 
Phys. Lett. B584, 48 (2004). 

A. Volya and V. Zelevinsky, Phys. Rev. C 67, 054322 
(2003). 

T. Berggren, Nucl. Phys. A 109, 265 (1968). 



[12] T. Berggren, Nucl. Phys. A 169, 353 (1971). 

[13] T. Berggren, Phys. Lett. B73, 389 (1978). 

[14] T. Berggren, Phys. Lett. B373, 1 (1996). 

[15] P. Lind, Phys. Rev. C 47, 1903 (1993). 

[16] G. Hagen, J. S. Vaagen, and M. Hjorth- Jensen, J. Phys. 

A: Math. Gen. 37, 8991 (2004). 
[17] B. Jonson, Phys. Rep. 389, 1 (2004). 
[18] N. Moiseyev, Phys. Rep. 302, 211 (1998). 
[19] R. R. Whitehead, A. Watt, B. J. Cole, and I. Morrison, 

Adv. Nucl. Phys. 9, 123 (1977). 
[20] K. Suzuki and S. Y. Lee, Progr. Theor. Phys. 64, 2091 

(1980). 

[21] K. Suzuki, Prog. Theor. Phys. 68, 246 (1982); 

[22] K. Suzuki and R. Okamoto, Prog. Theor. Phys. 92, 1045 

(1994); ibid. 93, 905 (1995). 
[23] S. Fujii, E. Epelbaum, H. Kamada, R. Okamoto, 

K. Suzuki, and W. Cockle, Phys. Rev. C 70, 024003 

(2004). 

[24] C. Buth, R. Santra, and L. S. Cederbaum, Phys. Rev. A 

69, (2004) 032505. 
[25] F. Chen, E. R. Davidson, and S. Iwata, Int. J. Quantum 



18 



Chem. 86, (2002) 256. 
[26] R. Santra and L. S. Cederbaum, Phys. Rep. 368, 1 
(2002). 

[27] S. Sack, L. C. Biedenharn, and G. Breit, Phys. Rev. 93, 
321 (1954). 

[28] P. Navratil and B. R. Barrett, Phys. Rev. C 57, 562 
(1998). 

[29] P. Navratil, J. P. Vary, and B. R. Barrett, 

Phys. Rev. Lett. 84, 5728 (2000). 
[30] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. C 

62 , 054311 (2000). 



[31] P. Navratil, G. P. Kamuntavicius, and B. R. Barrett, 
Phys. Rev. C 61 , 044001 (2000). 

[32] N. J. Higham, Num. Algorithms 15, (1997) 227. 

[33] E. D. Denman and A. N. Beavers, Appl. Math. Corn- 
put. 2, (1976) 63. 

[34] D. J. Dean and M. Hjorth- Jensen, Phys. Rev. C69, 
054320 (2004). 

[35] K. Kowalski, D. J. Dean, M. Hjorth- Jensen, T. Papen- 
brock, and P. Piecuch, Phys. Rev. Lett. 92, 132501 
(2004). 



